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ABSTRACT 



The hierarchical clustering inherent in A-CDM cosmology seems to produce many of 
the observed characteristics of large-scale structure. But some glaring problems still 
remain, including the over-prediction (by a factor 10) of the number of dwarf galax- 
ies within the virialized population of the local group. Several secondary effects have 
already been proposed to resolve this problem. It is still not clear, however, whether 
the principal solution rests with astrophysical processes, such as early feedback from 
supernovae, or possibly with as yet undetermined properties of the dark matter itself. 
In this paper, we carry out a detailed calculation of the dwarf halo evolution incor- 
porating the effects of a hypothesized dark-matter decay, D —> D' + I, where D is 
the unstable particle, D' is the more massive daughter particle and I is the other, 
lighter (or possibly massless) daughter particle. This process preferentially heats the 
smaller haloes, expanding them during their evolution and reducing their present-day 
circular velocity. We find that this mechanism can account very well for the factor 4 
deficit in the observed number of systems with velocity 10-20 km s _1 compared to 
those predicted by the numerical simulations, if Am/m^ ~ 5 - 7 x 10~ 5 , where Am 
is the mass difference between the initial and final states. The corresponding lifetime 
r cannot be longer than ~ 30 Gyr, but may be as short as just a few Gyr. 
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2 Majd Abdelqader and Fulvio Melia 
1 INTRODUCTION 

Observations of the cosmic microwave background (CMB) radiation with the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) have facilitated the precision measurement of several cos- 
mological parameters (Bennett et al. 2003, Spergel et al. 2003), including the mass-energy density 
of the Universe, Q,, which appears to be close (if not equal) to its critical value. Baryons contribute 
only about 4% of this; the rest is presumably in the form of dark matter (DM; roughly 22%) and 
dark energy (~74%). In addition, WMAP's detection of early reionization also rules out the pres- 
ence of a warm DM, so the non-baryonic component must be cold (CDM). Together with earlier 
observations by other finer scale CMB experiments, and with the Hubble Key Project (Mould et 
al. 2000), which provided the unprecedentedly accurate value H = 71 + 6 km s _1 of the Hubble 
constant, this combined body of work has lead to a consensus that the Universe is adequately de- 
scribed by the so-called flat A-CDM standard model, in which dark energy is the manifestation of 
a cosmological constant A. 

The existence of dark energy has been confirmed through the analysis of Type la supernova 
data (Riess et al. 1998; Perlmutter et al. 1999), providing even stronger evidence that this com- 
ponent of Q. has negative pressure leading to an acceleration of the Universe in the current epoc. 
To be fair, however, it is not yet entirely clear whether these results require a cosmology with a 
true cosmological constant, in which matter and radiation become dominant looking back towards 
redshifts z >2-3, or whether the current acceleration might be due to a so-called scaling solution, 
in which the dark energy density scales with matter, and affects the formation of structure even at 
early times (see Melia 2007, 2008, and references cited therein). 

Very little is known about dark matter, and almost nothing is understood about dark energy. 
Their nature is one of the biggest mysteries in contemporary physics. Yet their influence is evi- 
dently quite important in the formation of large scale structure (LSS). The hierarchical clustering 
inherent in A-CDM seems to produce many of the characteristics observed in the local (e.g., Murali 
et al. 2002; Abadi et al. 2003) and high-redshift (e.g., Springel, Frenk, and White 2006) Universe. 
In hierarchical models, smaller dark matter haloes on average collapse earlier than larger ones, 
when the density of the universe was higher (e.g., Kravtsov et al. 1998). The current mass func- 
tion, however, is determined not only by the halo formation history, but also by their merger rate 
which, over time, tends to deplete the dwarf-galaxy end of the distribution. 

But some glaring problems still remain. Numerical simulations involving collisionless CDM 
predict dark haloes with steep cusps in their centre (Navarro et al. 1996), whereas most of the 
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Decaying Dark Matter 3 

observed rotation curves of dwarf galaxies and low surface brightness galaxies indicate constant 
density cores (see de Blok et al. 2003, and references cited therein). A related problem is the over- 
prediction of the number of dwarf galaxies within the virialized population of the local group. 
CDM simulations over-predict the number of satellite galaxies orbiting a Milky Way-sized galaxy 
by a factor of 10 (Klypin et al. 1999, Moore et al. 1999, Diemand et al. 2007, and Simon & Geha 
2007). Generally speaking, both of these problems may be described as a A-CDM prediction of 
too much power on small scales. 

Since in hierarchical models smaller galaxies merge to make larger ones, the number of re- 
maining dwarf haloes is an important diagnostic to test both the hierarchical picture and the pro- 
cess of halo condensation in the evolving universe. Thus, the dwarf-halo deficit may be taken as 
an indicator of new physics associated with dark matter (and/or dark energy). For example, the 
small-scale power can be reduced by substituting warm dark matter for CDM. But as we have 
seen, WMAP observations have already ruled this possibility out. 

Other simple attempts to fix the dwarf-halo deficit problem are not well motivated physically. 
Some involve altering the fundamental nature of dark matter by introducing self-interaction, or an- 
nihilation (see, e.g., Spergel and Steinhardt 2000; Kaplinghat, Knox, and Turner 2000; and Giraud 
2001, among others). Without a proper indication from physical considerations, all of these mod- 
els contain free parameters tunable to fit the observations. For example, the self-interacting dark 
matter scenario rests on the viability of a huge velocity-dependent cross section. Unfortunately, 
the implied large interaction rate is inconsistent with most weakly interacting, massive particle 
and axion theories (see, e.g., Hennawi and Ostriker 2002). 

Of course, the dwarf-halo deficit may be due to reasons other than DM physics. Many of 
them may be invisible because they contain a very small amount of luminous matter, either be- 
cause of early feedback from supernovae (Dekel and Silk 1986; Mac Low and Ferrara 1999), or 
because their baryonic gas was heated by the intergalactic ionizing background radiation (Rees 
1986; Barkana and Loeb 1999). Others may have turned into high-velocity clouds in the Local 
Group (see, e.g., Blitz et al. 1999). The deficit may not even be real if the Universe is actually 
older than its current inferred age, which in reality is only the light-travel time to the cosmic hori- 
zon rather than the Big Bang (Melia 2007, 2008). In such a scenario, the dwarf haloes would have 
had more time to merge, depleting the lower end of the mass function and bringing it into better 
alignment with the observations. 

In any case, it is still too early to tell whether or not the discrepancy between the cosmological 
simulations and observations really indicates a major problem for hierarchical models in A-CDM. 
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4 Majd Abdelqader and Fulvio Melia 

Several of the effects we have listed here may resolve at least part of the deficit problem. However, 
given that this is still an open question, observations of the halo mass function also allow us to 
explore non-astrophysical reasons for the discrepancy, with the goal of learning more about the 
nature of dark matter. 

Our focus in this paper is the suggestion that DM particles may be unstable to decay (see, e.g., 
Davis et al. 1981; Turner et al. 1984; Turner 1985; Dicus and Teplitz 1986; Dekel and Piran 1987; 
Sciama 1990; Cen 2001a; Sanchez-Salcedo 2003). The impact of the interactions we describe 
above and/or decays is almost always to provide a mass-dependent expansion of the cusps and 
haloes to lower the core density and to reduce the number of small galaxies. However, attempts 
to couple these ideas to particle physics have been few and ambiguous, partly because these have 
been concerned more with global effects, rather than aimed at finding specific particle properties 
that may be consistent with the requirements to fix the problem. Our goal in this paper is to begin 
a more careful search for the properties DM particles must have in order to account for the deficit 
of dwarf haloes, if in the end their decay is indeed responsible for the observed effect. 

Our approach here is closest in spirit to the work of Cen (2001a) and Sanchez-Salcedo (2003), 
though their papers had different goals. Cen's (2001a) primary interest was to demonstrate that if 
DM particles decay and become relativistic, they escape the virialized halo, whose remaining en- 
ergy then exceeds that required to sustain virial equilibrium and forces it to expand. His suggestion 
was that the overproduction of dwarf haloes may be solved not by removing them, but by modify- 
ing them into failed, dark galaxies, in which star formation has been quenched due to the effects 
of evaporation and expansion. This is an intriguing idea, though not yet sufficiently developed to 
provide a useful probe into the properties of the particles themselves. 

Sanchez- Salcedo's (2003) goal was to mitigate the disparity between the very steep central 
cusps in dark haloes of dwarf galaxies predicted by A-CDM and the relatively flat distributions 
actually seen in these systems. He demonstrated that if DM decays into a relativistic, nonradiative 
light particle, plus a stable massive particle with some recoil velocity in the center-of-mass frame, 
the former escapes the bound system while the latter remains with an energy exceeding that of the 
parent, forcing the halo to expand. 

In this paper, we introduce several new aspects of the DM-decay scenario, including the time- 
dependent and mass-dependent halo formation and hierarchical-merger rates in order to more ac- 
curately gauge the impact of heating on the present-day circular velocity distribution. There are 
too many aspects of the DM decay to consider in just one set of calculations, so we here restrict our 
attention to cases in which at least one of the decay products remains within the halo, maintaining 
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its mass, though heating it with the liberated energy. Other regions of the DM particle phase space 
will be explored elsewhere. 

In the next section, we summarize the circular-velocity data and demonstrate the nature of the 
dwarf-halo deficit problem. Then in §3 we describe a technique for following the mass-dependent 
formation and destruction of haloes as the universe evolves. In §4 we describe the DM-particle 
decays and how we incorporate the impact of this process into our calculational algorithm. We 
present the results of our calculations in §5, and discuss them in §6. 

2 THE OBSERVED CIRCULAR- VELOCITY DISTRIBUTION OF DWARF HALOES 

The dwarf-galaxy deficit was first quantified when the observed number of dwarf galaxies in the 
local group was compared to high-resolution cosmological simulations of Klypin et al. (1999) 
and Moore et al. (1999), under the reasonable assumption that the local group is an adequate 
representation of what is happening throughout the cosmos. If we assume that each small dark 
matter halo contains a dwarf galaxy, then there is a substantial discrepancy between theory and 
observation. At the time of these simulations, there were only 13 known satellites of The Galaxy, 
while both simulations predicted roughly 10 times that number of satellites for a Milky Way-sized 
halo (Klypin et al. 1999, and Moore et al. 1999). 

In the last several years, the number of pertinent observations has increased substantially, and 
new cosmological simulations have been completed with substantially higher resolution than the 
original calculations. The Sloan Digital Sky Survey (SDSS hereafter) has uncovered 10 more 
Milky Way companions (Belokurov et al. 2007). The discrepancy between the latest observations 
and the most recent simulations still exists, however, and it appears to be a function of mass. In 
the range of primary relevance to this paper (i.e., 10-20 km s _1 ), the disparity is a factor 4 when 
compared to the latest high resolution N-body simulation known as Via Lactea (Diemand et al. 
2007), even after weighting the new dwarf galaxies by a factor of 5 to account for the limited 
coverage of the SDSS; at 6 km s _1 , the discrepancy increases to a factor 10 (Simon & Geha 2007). 

3 FORMATION OF BOUND OBJECTS IN THE HIERARCHICAL CLUSTERING 
SCENARIO 

For reasons that will soon become apparent, the impact of decaying dark matter on the evolution of 
haloes depends on their formation history. For simplicity, we here use a semi-analytical approach 
that describes the halo formation rate as a function of mass and time, optimized to reproduce 
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6 Majd Abdelqader and Fulvio Melia 

numerical simulations of structure formation. A good starting point is the Press-Schechter mass 
function (hereafter PS), or one of its modified forms, which is a reasonable representation of the 
overall halo distribution resulting from these numerical simulations (Press and Schechter 1974). 
However, PS is a number density that combines both the formation and merger histories of the 
haloes, so it does not provide their formation rate explicitly. None the less, a formation rate can 
be extracted from the PS formalism by taking the comoving time derivative of the mass function, 
identifying the term corresponding to the formation rate, and multiplying it by the survivability 
probability. This procedure subtracts the halo destruction rate, and is necessary for our purpose 
since we only want to consider haloes that formed in the past and survived to the present without 
merging with others to form even bigger structures (see, e.g., Sasaki 1994; Kitayama & Suto 1996). 
The formation rate can be written as (Sasaki 1994) 

1 da(t) 6 2 (t) 
N form (M, t) dM dt = ——LLn ps (M, t)-^- dM dt , (1) 
a(t) dt cr 2 (M) 

where M is the mass of the formed gravitational structure, t is the formation time measured in 
comoving coordinates starting from zero at the Big Bang, a(t) is the cosmological scale factor 
normalized to unity at the present epoch t , N PS (M, t) is the Press-Schechter mass function, 8 c (t) = 
^ is the critical density threshold for a spherical perturbation to collapse by time t (6 C - 1 .69 for 
Q.0 = 1), and cr(M) is the rms density fluctuation smoothed over a region of mass M. However, 
we are interested in haloes surviving to the present epoch, thus the above formation rate function 
must be multiplied by the probability p(t\, t 2 ) that an object which exists at t\ remains at t 2 without 
merging, which is given by p(h,t 2 ) = a(ti)/a(t 2 ) (Sasaki 1994). The formation rate distribution 
of surviving haloes at the present epoch becomes F(M, t) dM dt = Nf orm (M, t) x p(t, t ) dM dt, 
which can be written explicitly as 

F(M, t) dM dt = A^ — M^ 12 exp 
dt a(ty 

where A is a normalization constant, n is the power-law spectral index chosen to be -2.5, following 
Klypin et al. (1999), and M Ci0 is the present mass scale of the knee taken to be 10 15 M©. 

In this paper, we adopt a flat universe Q = 1> represented roughly as Q m>0 = 0.3, and Q Aj0 = 
0.7, which gives a(t) = (0.3/0.7) 1/3 sinh 2 ' 3 (1.21 [t/tO]), normalized to unity at the present epoch 
to = 13.7 Gyr. Figure 1 shows the formation rate as a function of time for two illustrative halo 
masses. The larger mass always has the smaller formation rate, which also peaks at later times, 
though this difference is not very obvious from the figure as the two are relatively close. 
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Figure 1. The relative halo formation distribution as a function of cosmic time for haloes surviving to the present epoch, for two illustrative masses: 
M = 1 x 10 8 M o and M = 3 x 10 8 M S . 

4 DM PARTICLE DECAYS AND THEIR IMPACT ON DWARF HALO EVOLUTION 

In this paper, we focus our attention on dark matter decay scenarios in which the parent particle, 
denoted by D, decays into one stable massive particle (£>') with a mass close to that of its parent, 
and one very light (possibly even massless) particle (€), following the notation of Sanchez-Salcedo 
(2003). Due to the extreme mass ratio between the daughter particles, the light particle carries 
most of the energy released by the decay and becomes relativistic, escaping the halo. The massive 
daughter particle remaining in the halo has an (average) energy very close to, but slightly larger, 
than that of its parent, forcing the halo to expand adiabatically. Since the kinetic energy of the light 
particle is much larger than its rest mass, we neglect the latter. To this level of approximation, the 
total mass of the halo remains unaffected by the decay because m D , w m D . 

The decay is represented as D — » D' + t. In the center-of-mass frame (denoted by a * super- 
script), i.e., in the rest frame of the parent particle D, the initial four-momentum vector before the 
decay is 

p*=(m D c,6), (3) 
and after the decay the final four-momentum is 

P} = P* D ,+P* C =( Y D m D ,c , -y* D m D ,v* Di , , ) + ( E]lc , E*/c ,0,0), (4) 

where we have put the direction of p c along +x, and along —Sc. 
Conservation of momentum leads to the equation 

E* = (m 2 D -ml)c 2 /2m D . (5) 

To find the particle's energy E c in the halo's rest frame, we need to Lorentz boost the physical 
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8 Majd Abdelqader and Fulvio Melia 

quantities using the velocity v = v D and the angle 0* between v and p r This gives 

E { = y D E*(l+/3 D cos9*). 
Averaging over all solid angles, we find that 

<E,> = -L J E e dQ = y D E*. 



We next define the unities s parameter 



X = 



Am 



m 



(6) 



(7) 



(8) 



the ratio between the change in rest mass Am = m D - m D , and the mass of the heavier daughter 
particle D' . The energy of the lighter (relativistic) particle, averaged over all angles 9*, becomes 



(E t ) = y D m D 



X(2+X) 



(9) 



'2(1 +x) 2 ' 

If both particles remain in the halo, then the change in the halo's energy for each decay is 
simply Am c 2 . However, the lighter particle escapes, so its energy must be subtracted. We find that 
for this type of decay, the rate of change in the halo's energy is therefore 

'■'halo . 2 /F \ l[ X X( 2 +X)\ 



dEi, 



dN 



(10) 



where N is the number of unstable particles. But the decay rate is 

<£L = ± NQe '^lr = _^ g -(»«/)/r t (11) 

dt dt t 
where tf is the time at which the halo forms, t is the time elapsed since the formation of the halo, 

r is the mean lifetime of the parent particle D, and No is the initial number of parent particles at 

the time the universe began its expansion. Combining these quantities, we can now find the rate at 

which the halo's energy changes with time: 



cLE 
dt 



dE_ 

In 



dN 



dt 



= (Nom D )c 2 (t+tf)/T I X X(2+X) 
r ^ \l+X ° 2(1+^)2 



(M)c 2 



-(t+tf)/T 



X 



7 D : 



X(2+X) 



(12) 



1+X " J 2(\+ X ) 2 > 

By knowing the rate of energy change, we can in principle find the change in size of the halo 
by expanding it adiabatically. However, in order to do that, we need to know its initial density 
profile. At the time of formation, we will assume the halo has a Navarro, Frenk, & White density 
profile (1997, NFW hereafter) 



PnfwO) 



Pc„, (r/r s )(l + r/r s ) 2 ' 
where p is the critical density, r s is the scale radius of the NFW profile, and 



(13) 



200 



6r] 3 In(l+i7)-i7/(l+i7) 



(14) 

© 0000 RAS, MNRAS 000, 000-000 



Decaying Dark Matter 9 

is a characteristic (dimensionless) density in terms of rj = r v /r s (the concentration parameter) and 
the virial radius r v (defined as the radius of a sphere of mean interior density 200p crit ). Although 
the NFW profile is an analytic function, it cannot be readily incorporated into our semi-analytical 
model. First, its distribution function f(r, v) cannot be obtained analytically, and must be found 
numerically. Furthermore, to find how the scale radius r s evolves with time, we need to find the 
halo energy as a function of r s , but this is not easy to do with the NFW profile. Assuming the halo 
is virialized, E tot = (K) + (U) = -\ (U) + (U) = \{U). The total potential energy is given as 
(Binney & Tremaine 1987) 

U = \J 3>(x)p(x)J 3 x, (15) 

where O is the gravitational potential. 

To evaluate this integral, we need to find <D(x) from p(x), which we can do by first finding the 
mass 

r 

M(r) = j Anr' 2 p(r') dr' . (16) 

o 

For the NFW profile, 

r/r s 



M NFW (r) = M -P^ , (17) 

NFwV ln(l + jy) - JL ' v 



_ln(l+r/r,) ^ 
n 

where M is the virial mass of the halo contained inside the virial radius. By definition, VO(x) = 
-F(x), and for the simple isotropic case, F(r) = -GM(r)/r 2 , which leads to 

* (r) = J5W) d/ . (18) 

n 

For an NFW halo, this give^ 



r ln(l + ?]) - j± 



Thus, solving for U in equation (|T5I) . we find that 



(20) 



GM 2 

'^7 



l - l/(l + rf) 2 - 2ln(l + ?/)/(l + 77) 



[ln(l +77) -77/d +77)f 

In this equation, the concentration parameter 77 changes with r s , which makes it difficult to find an 
analytic expression for drjdt in terms of dE/dt. 

For these reasons, it is convenient to translate the NFW profile into an equivalent Plummer 



To get the correct energy, the potential <D(x) must be adjusted with an additive constant such that <£(f", ) = —GM/r v since the NFW halo density 
must vanish for r > r v . 
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distribution, 



which is mathematically easier to evolve in time. In this expression, r p is the Plummer scale radius. 
A principal virtue of the Plummer profile is that it solves the Lane-Emden equation for a self- 
gravitating, polytropic gas sphere. The corresponding distribution function may be written 



f P (r,v) = B 



4 



GM 1 , 
— - -v 



r 2 + r 2 „ ^ 



7/2 



(22) 



where B is a normalization constant, and it is trivial to show from equations (fl6l) and (fT8l) that 



r 3 



^ = M ( 2 ^ ' (23) 
and 

-GM 

<W = • (24) 

V r p + r2 

Thus, the total energy of a Plummer halo (from equation IT3T) is 

64r p 

In making the transition from an NFW profile to its corresponding Plummer form, we use the 
physically reasonable criterion that two haloes should have the same mass and energy. Therefore, 
equating equations (1211 and (|25T) . we get 

[ln(l + 77) - 77/(1 + 77)] 



3nf r — ■ - ■ -* 2 



^ = r ' 16 



(26) 



l-l/(l+77) 2 -21n(l+77)/(l+77)J 
Following Navarro, Frenk, and White (1997), we will further assume that rj = 20 at the time of 
formation (for the smallest haloes), though our results are insensitive to its actual value. 



Figure 2 illustrates the differences in circular velocity (v C i rc = y/GM(r)/r) for an NFW halo 
(with mass 3 x 10 s M Q ) and the corresponding equivalent Plummer form. The Plummer halo may 
not fit the observed velocities as well as NFW, but their maximum circular velocities are almost 
the same, and since the Plummer model permits us to obtain an analytic solution for the halo's 
evolution in time, we will use it in all our calculations under the assumption that the collective 
behavior of self-gravitating virialized haloes is similar for slightly different profiles. Nevertheless, 
we will still need to use the NFW profile to obtain the initial characteristics of the halo at the time 
of its formation. 

Now, using the Plummer distribution (equation [22]), we can simplify equation (fT2l) by evaluat- 
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Figure 2. Circular velocity curves of two haloes with the same mass (M = 3 X 10 Mo) and the same total internal energy, but one with a Plummer 
density profile, and the other with an NFW density profile. Radii are in units of the initial Plummer scale radius, r„ j, and circular velocities in units 
of kilometers per second. 



ing the average rate 

( 



'?)-/§«--"'" 

_(M)c 2 c _ (t+tf)/T \ X 



l+X I 1 + 64c 2 rJ 2(1 + X ) 2 



3nGM\ X (2 + x) 



And combining equations (1251) and (1271) . we therefore get 



dr p 64 rj dE 



dt 



37TGM 2 dt 

r 2 

P 1 „-(t+tf)/T 



64 rl c 2 



3nGM t 



X 



3nGM\x(2+x) 



1+* 



/ 3nGM\ 
\ + 64c%j 



64c 2 r p ) 2(1 +x) 2 



This differential equation can be solved analytically, yielding the result 



r P (t) 



(27) 



(28) 



• p.i 



64 X c 2 r pJ exp [-^^//-] + (3nGM(x + 2) - 64^c 2 r p ,-) exp [-^ e -(«/)/r] " 



(29) 

where r p j is the initial Plummer scale radius at the time of formation. Note that the impact of DM 
decay on the evolution of the halo does not depend on the mass of the individual particles D and 
D', but rather on the ratio of these masses represented by the unities s parameter^, as well as on 
the mean lifetime r of the parent particle D. There is an underlying assumption here that the halo 



2 Note that y D = 1/ y" 1 - v 2 /c 2 a 1 + v 2 /2c 2 since particle D is non-relativistic. 
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Figure 3. Circular velocity curves at different times (indicated) of an expanding halo formed at f = \Gyr with mass 3 X 10 8 M o . The expansion is 
due to the decay of CDM particles with^ = 4x 10~ 5 and r = 5 Gyr. The radius is in units of the initial Plummer scale radius, and circular velocities 
are km s . 



keeps a Plummer profile throughout the expansion. When comparing the decay time scale, defined 
as t decay = \E/dE/dt\, to the dynamical time t dyn = -\j3nJT6Gp (Binney & Tremaine 1987), we find 
that tdecayltdyn > 10, implying that the halo always equilibrates to its virialized profile fast enough 
to justify the quasi-equilibrium approximation. Figure 3 illustrates the expansion of a halo formed 
at tf = 1 Gyr, showing how the maximum circular velocity decreases with time. 



5 COMPUTATIONS AND RESULTS 

To explore the global impact of our model on the halo distribution, we incorporate the effects of 
the DM decay on each individual halo (equation [29b and its formation rate (equation [2]) through 
a series of calculations. We assess the consequences of this process by examining the modifica- 
tions to the distribution of maximum circular velocity (v cl > c hereafter) under two assumptions: (i) 
that at the time of formation, all haloes with the same mass will have the same concentration pa- 
rameter. Although this is evidently incorrect for individual haloes, we are considering the global 
behavior, for which an average concentration parameter will suffice. Thus, in the case of no decay, 
all virialized haloes with the same mass will have the same circular velocity regardless of when 
they formed, since their concentration parameter does not change with time and only depends on 
the halo's mass; (ii) that the stellar dispersion velocity is proportional to the maximum circular 
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Figure 4. Circular velocity distribution of DM haloes for a fixed mean lifetime r = 5 Gyr, and different values of the decay parameter % = Am/m D i . 
The circular velocity distribution without decay is provided for comparison. 

velocity of the host dwarf haloid Thus, with DM decay, the haloes expand and their concentration 
will correspondingly change with time, which in turn alters the circular velocity and the observed 
stellar dispersion velocity. So haloes with different masses that formed at different times in the 
past may end up with the same circular velocity in the current epoch. 

For a given decay parameter ^ and a mean lifetime r, the evolution of the halo still depends on 
when it formed and on its mass. For example, if the halo formed late relative to the mean lifetime 
t, then most of the unstable DM particles will have already decayed by then, and the halo will 
therefore experience no significant expansion. Furthermore, according to equations (fl"2)) and (T2~5l) . 
dE/dt oc M, while E oc M 2 , which means that (dE/dt)/ E oc 1/M. This is a crucial dependence 
of this process on mass since it guarantees a relatively stronger impact on the smaller haloes. For 
example, if a small halo with M « 10 8 M Q experiences a significant expansion for a given set of 
decay parameters (for which, say, (dE/dt)/ E is of order unity), then a much larger halo, e.g., with 
M~ 10 10 M o , will be unaffected by the decay since now (dE/dt)/ E « 10~ 2 . 

Figure 4 demonstrates the impact of changing the decay parameter x on the circular- velocity 
distribution in the range 8-22 km s _1 , given a fixed mean lifetime. For this calculation, the forma- 

3 This assumption has been extensively used to infer the circular velocity for elliptical galaxies and dwarf spheroidals from the observed stellar 
dispersion velocity v clrc = V3cr, where a is the observed stellar dispersion velocity (Klypin et al. 1999, Simon & Geha 2007) 

© 0000 RAS, MNRAS 000, 000-000 



1 4 Majd Abdelqader and Fulvio Melia 




Figure 5. Circular velocity distribution of DM haloes for a fixed decay parameter ^ = 4x10 5 , and different values of the mean lifetime r. The 
circular velocity distribution without decay is provided for comparison. 

tion time distribution function was normalized such that without any DM decay, 100 haloes would 
be produced with v c \ rc between 10 and 20 km s" 1 . We can see that the decay decreases the number 
of haloes with lower circular velocities more than those with higher ones. Furthermore, the impact 
of the DM decay on the velocity distribution increases with^, which is not surprising given that 
a higher x corresponds to a higher recoil velocity of the remaining particle D' and, therefore, a 
greater expansion. It is important to note that our approximations cease to be valid at^ « 7 x 10~ 5 , 
because for higher values the recoil velocity of D' becomes comparable to its escape velocity in 
the haloes we are considering. For such high values of x, some of the D' particles would start 
escaping right after the decay, depending on the velocity and initial position of the parent particle 
and on the angle 9* that the recoil velocity makes with the velocity of the parent particle. As a 
result, the mass of the halo would decrease with time, an effect that is not being taken into account 
right now. 

We have also examined the impact of varying the mean lifetime r on our model, the results of 
which are summarized in figure 5, for a fixed value of x- For relatively small values of r (« 0.1 
Gyr), the decay has virtually no impact on the circular velocity distribution since most of the DM 
particles will have already decayed before the vast majority of the haloes formed. In these cases, 
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Figure 6. Contour plot representing the effectiveness of the DM decay in reducing the number of dwarf haloes with v c irc between 10 and 20 km 
s , given a set of decay parameters x - Am/m D > and r. The shaded regions represent the ratio of the number of haloes predicted without decay 
(for v c irc between 10 and 20 km s ) to that obtained with DM decay. 



the haloes experience no significant expansion. As r increases, however, most of the DM particles 
decay after the majority of the haloes have formedO The impact of DM decay reaches a maximum 
"effectiveness" in decreasing the number of haloes for r ~ 10 Gyr. As one would expect, larger 
values of r produce less significant results since most of the unstable particles would not have 
decayed by the present time. 

The full dependence of the circular-velocity distribution on DM decay is shown in figure 6, 
which includes the effects of both^- and r. The contour levels represent the effectiveness of the 
decay in decreasing the number of dwarf haloes within the velocity range 10-20 km s _1 . The 
various shaded regions represent the expected number of haloes without decay divided by the 
corresponding number when expansion occurs with a given set of parameters x an d t. 

Besides adjusting the predicted circular velocity distribution to bring it in line with observa- 
tions, the most important observational signature of DM decay is the dependence of the concen- 
tration parameter of dwarf galaxies on redshift. At relatively large redshifts, the haloes would have 
formed recently and a smaller fraction of the DM particles would have decayed, so neither the 
concentration parameter, nor the stellar dispersion velocity, would have been influenced greatly. 
This is to be contrasted with what would have happened to haloes observed in the current epoch. 
Therefore, distant dwarf haloes at large redshifts would be expected to be brighter (on average) 



4 This implies that the dependence of our model on r is strongly related to the formation time distribution function. Thus, using a different formation 
function can produce a qualitatively different dependence on r. 
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and to have higher dispersion velocities compared to their nearby (lower redshift) counterparts 
with similar masses. 

6 CONCLUSION 

We have shown that the decay of unstable CDM particles can fully account for the deficit of 
dwarf galaxies in the local group, and have identified some of the particle properties required to 
achieve this result. In figure 6, the lifetime r and energy conversion fraction x must have values 
consistent with the thick black line between the two gray regions in order to reduce the number of 
dwarf haloes within the range of velocities 10-20 km s _1 by a factor ~ 4, in agreement with the 
observations (Simon & Geha 2007). Broadly speaking, the DM decay model works very well as 
long as x ~ 5 - 7 x 10 5 . The lifetime cannot be longer than ~ 30 Gyr, and may be as short as 
only a few Gyr, depending on the precise value of x- Note, however, that we have here restricted 
our analysis to cases in which the more massive decay particle D' remains bound to the halo. We 
may find other regions of x ~ t phase space that produce reasonable results when this restriction 
is removed. 

We emphasize that this model can reduce power on small scales consistent with the observa- 
tions without altering the number of Milky Way-sized galaxies. Very importantly, we have shown 
that although the expansion produced by these decays changes the circular velocity distribution, it 
does not change the halo mass function, at least not directly. It is beyond the scope of the present 
paper to seek the ultimate fate of dwarf haloes expanding to circular velocities below 10 km s _1 , 
which remains an open question, though several possibilities have been proposed by Cen (2001a, 
2001b). 

It should be pointed out, however, that in addition to expanding (preferentially) the smaller 
haloes, and thereby reducing their measurable velocity dispersion, DM decay would also have the 
effect of speeding up their rate of evaporation within the Milky Way's tidal field. In this way, dwarf 
haloes would be removed from the overall velocity distribution, not only due to their migration in 
velocity space towards the low end, but would also be removed entirely due to evaporation. This 
effect would not be evident with the larger haloes, for which DM decay would have little impact 
on their condensation (and hence on their velocity profiles). 

Interestingly, some limits on DM-decay models have already been established based on limits 
provided by the cosmic y-ray background. Although it is beyond the scope of the present paper to 
consider the implications of our work on all possible DM scenarios, it is useful to see how coupling 

© 0000 RAS, MNRAS 000, 000-000 



Decaying Dark Matter 17 

our astrophysically-motivated simulations to the various particle physics proposals could develop 
in the future. For example, in their consideration of WIMPs decaying to Kaluza-Klein gravitons 
and gravitinos, Cembranos et al. (2007) demonstrated that both the energy spectrum and flux of the 
observed diffuse MeV y-ray excess may be explained by decaying DM with ~MeV mass splittings. 
In this picture, a decay timescale of 10 Gyr would require a mass splitting of ~ 10 MeV, for which 
the unstable DM particle would then have a mass ~ 1 TeV within the context of our model. In a 
second example, Kribs and Rothstein (1997) placed bounds on long-lived primordial relics using 
measurements of the diffuse y-ray spectrum from EGRET and COMPTEL. They concluded that 
relics decaying predominantly through radiative channels are excluded for lifetimes between 10 5 
and 10 15 years. Since the DM decay timescale in our model fits within this range, the radiative 
decay of relics such as these could not resolve the dwarf spheroidal problem. 

We have kept our analysis semi-analytical in order to better gauge the impact of our assump- 
tions and chosen parameter values. Of course, to get a more accurate set of results, one should 
couple the properties of decaying DM particles with a more realistic N-body simulation. We in- 
tend to carry out such a calculation in the near future and will report the results elsewhere. 
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